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Abstract 

It is shown that in a Minkowski space of total space-time dimension 
D = d + 1, the orbits of the planetary motion are stable only if the 
total dimension of space-time is D < 4. The proof is performed in a 
fully didactic way. 



1 INTRODUCTION 

Bertrand's theorem establishes that the only central potentials with stability 
are the harmonic oscillator and the Coulomb's one [2]. It has been related 
the closeness of the classical orbits to their symmetries [3]. The existence of 
different kinds of factorization operators for the radial Schrodinger equation 
corresponding to the Kepler and harmonic oscillator problems has been at- 
tributed to the classical stability of these systems according to the Bertrand's 
theorem [4]. Also, a generalization of this result has been developed when 
the potential has suitable extra angular momentum [5]. More subtle points 
has been studied like the connection between closed orbits and their stable 
symmetries for small perturbations [6]. It has been found that when all eight 
symmetries are stable the orbits are closed, but not necessarily if only three 
symmetries of them are stable, as it happens in the relativistic case where 
the perihelion of the orbit advances [6]. 
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Recently it has been studied the circular orbits stability for the Schwarzs- 
child geometry (D = 4) in commutative and noncommutative spaces [7]. 
Schwarzschild solution in commutative space possesses stable circular orbits 
for a radial distance r > 6GM and unstable ones for 3GM < r < 6GM, 
where M is the mass of the attractor center and G is the gravitational con- 
stant. The space non-commutativity increases the radius of stable circular 
orbits [7] . Gravitational classical instability in higher dimensions (D > 4) has 
been explored for various static and non static spacetimes: black holes, brane 
worlds, generalized Schwarzschild, AdS black strings, "bubbles of nothing" 
[8] and hyperbolic spaces [9]. 

The aim of the present paper is to show that the planetary circular orbits 
in a Minkowski spacetime of total dimension D = d+1, are stable only when 
D < 4. To obtain this result the Bertrand's theorem plays a central role 
which allow us to give an elegant and didactic proof. In section 2, we obtain 
the gravitational force in d spacial dimensions due to a particle of mass M. 
In section 3, we sketch the Bertrand's theorem, and it is applied to constraint 
the spatial dimensions in order to obtain the stability of circular orbits. In 
section 4, we give the concluding remarks. 



2 GRAVITATIONAL FORCE IN A 
^-DIMENSIONAL SPACE 

The magnitude of the attractive gravitational force between two point par- 
ticles m and M in the three-dimensional space is given by 

F=-G^ f , (1) 
which satisfies the Gauss' law (setting m — 1) in three dimensions 

/ F-dd=-4nG^M enc , (2) 

JS 2 (r) 

where M enc is the enclosed mass by the Gaussian sphere S 2 and G^ = G. 

We are interested in generalizing expression (1) to a d- dimensional space. 
In order to obtain this, we demand that the gravitational force satisfies the 
Gauss' law (2) for any d value, remaining fixed the value of the right hand 
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side. This means that the flux of the gravitational force on any sphere S d 1 (r) 
is always equal to the enclosed mass times 4nG^ [1]. 

Some remarks about da are necessary. If we consider the spheres S 1 
and S 2 , their areas are equal to 2nr and 4nr 2 , respectively. The area 
of the sphere S 3 is some alike a three-dimensional volume. To higher- 
dimensional spaces the area of a sphere S 4 ' 1 of radius r is known as the 
volume, a =(vol(S ,d ~ 1 (r))). In what follows we use this terminology to de- 
note areas for any d value. Hence, the volumes for S 1 and S 2 are 

vol^r)) = 27rr - ( 3 ) 
vol(,S 2 (r)) = 4vrr 2 . (4) 

Thus, the factor 4n in equation (2) is equal to the volume of the unit sphere 
S 2 . Equations (3) and (4), lead to infer that for a sphere in a rf-dimensional 
space the equality vol(S ,( - c!_1 - ) (r)) = r d ~ 1 vo\(S d ~ 1 ) must hold, being vol(<S' d-1 ) 
the volume of the unit sphere. Moreover, the volume differentials (r fixed) 
are 

d[vol(S ,1 (r))] = rd6 = rd[vo\(S% (5) 
d[vo\{S 2 (r))] = r 2 sin(0)d0# = r 2 d[vo\{S 2 )}, (6) 

from which we conclude that in general, for a <i-dimensional sphere, the 
equality 

da = dlvoliS*- 1 ^))] = r d - 1 d[vo\(S d - 1 )] (7) 

holds. 

In order to the gravitational force flux on the sphere S' d ~ 1 (r) be the fixed 
constant AttG^ 3 \ according to the result (7), the force between two massive 
point particles must depend on l/r d_1 . We propose it has the form 

F^-G^r. (8) 

Now we focus our attention on the right hand side of equation (2). The mass 
M as an intrinsic matter property must not depend on the space dimension. 
Hence, the product 4^(7^ = vo\(S 2 )G^ must have a fixed value for any d. 
Since in higher dimensions the area (volume) of the unit sphere depends on 
the space dimension, the gravitational constant must depend on the dimen- 
sion too. Hence, vol(S' c( ~ 1 )G'^- ) = AwG^. Above remarks allows to write the 
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generalized Gauss' law as 

f F-da = -vo\(S d - 1 )G^M enc . (9) 

which is satisfied by the gravitational force 

F — — 47rG(3) TE f (10) 

for m — 1. 

The next step is to find the volume of a unit sphere of dimension d — 1. 
We proceed as in [10] and calculate the integral 

/■oo /-oo ;-oo 2 f 2 

i ' ' ' / dx\dx2 ■ ■ ■ dxae~ r = / e~ r dt> (11) 
-oo J— oo J— oo Ji? d 

in a// space in two ways. 

First, since r 2 = + x\ + ... + a; 2 ,, the integrand of (11) is separable and 
it reduces to a product of one-dimensional integrals. Thus, 



However, each integral of the product is well known from the probability 

pOO T o>2 i ml r t A 

-oo' 

?d 



U dxte-t (12) 
i=i J "°° 

roduct is well kno 
theory, it is f^dxie - ** = . Therefore, Id = . 

Second, the value of Id is obtained by considering the space R d as divided 
by thin spherical shells. At fixed r the space is the sphere S' rf_1 (r) and the 
volume of the space between r and r + dr is dv =vol(S d " 1 (r))dr . This fact 
allows to write equation (11) as 

roo „ 

I d = / vo\{S d -\r))e- r2 dr (13) 

roo „ 

= volf^- 1 ) / r d - l e- r dr (14) 

JO 

= -volfS"*" 1 ) / t^e-'dt, (15) 
2 Jo 

where we have used the relation vol(S ,( - d_1 - ) (r)) = r d ~ 1 vo\(S d ~ 1 ) and per- 
formed the variable change t = r 2 . The last integral is identified with the in- 
tegral representation of the gamma function defined by T(x) = / °° t x-1 e - *dt, 
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with x > 0. With this result we obtain that I d = |vol(S' fl! " 1 )r(f ). Thus by 
equating our two results for Id, we obtain 



Results (10) and (16) complete our purpose to find the gravitational force 
between two point particles in a d- dimensional space. 

Bertrand's theorem ensures that the circular planetary orbits around the 
sun in our four- dimensional space-time are stable under small perturbations 
[2], and relates its stability to the values of the power in potentials with the 
form V = r u . We notice that the force (10) due to one particle of mass M 
on a particle of mass m is proportional to \jr d ~ x . This fact will allow us 
to apply Bertrand's theorem to find the values of the spatial dimension d in 
order to have stable planetary orbits in higher dimensions. 



3 THE BERTRAND'S THEOREM AND 
STABILITY OF THE CIRCULAR ORBITS 
IN A ^-DIMENSIONAL SPACE 



We begin by reviewing the Bertrand's theorem proof. Since the total angular 
momentum is conserved, this allows us to set an inertial frame in the mass 
center of the system. Also, since the force is central, the angular momentum 
is conserved and restrain the orbit to be plane. Hence, the motion of the 
two particles is reduced to that of a single one with reduced mass /x = ™+ M 
under the potential V(r), being r the distance between particles. It is usual 
to set the z-axis in the angular momentum direction, orthogonal to the orbit 
plane. Also, since V(r) does not depend on time t, the total energy E is 
conserved. 

Since the lagrangian of the system is 



vol^" 1 ) 




(16) 



L = T-V = -fi(r 2 + r 2 6 2 ) - V(r) 
2 



(17) 



and does not depend on the variable 6, the angular momentum £ = fir 
conserved. The total energy of the system is 



E = -/i(r 2 + r 2 6 2 ) + V(r) 



(18) 
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2^ + - 2 - 2 +V(r) 



,> 2 + V ef , 



(19) 
(20) 



from which is obvious the definition of the effective potential V e f(r). For the 
case in our study, the reduced mass is ji — to. Thus, it follows that 



dV f 



ef 



dr 



or 



where F = — 

or 



ef 



dV(r) 
dr 

= F + 



m+M 

mr 3 



(21) 
(22) 



The condition to have circular orbits is imposed by r = 0, and they are 
stable only when the potential has an effective minimum at a distance r = r . 

must be equal to 



sv. 



>-f 



r=r 



This fact implies that for circular orbits, F e f = — ih 

zero (whenever the case, a maximum or a minimum). Thus, equation (22) 
reduces to 



F = 



(23) 



By demanding that the effective potential V e f(r) get a minimum at r = r , 
from equation (21) it must satisfy 



d 2 V e 



ef 



this implies 



dr 2 

dF 

dr 



d 2 V 



r=ro 



dr 2 



+ 3- 



r=ro 



mr q 



> 



> 



r=r 



mrQ H 



-Fir) 
r 



\r=r i 



or 



dF 

dr 



r=r 



< F(r ). 

r 



(24) 



(25) 



(26) 



This inequality must be satisfied in order to have stable circular orbits at 
r = r . 
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From equation (10) and (16), the force between two particles in a d- 
dimensional space is given by 



~ "vol^- 1 )^- 1 ~ vrl-i r d - 1 ' ( ~' } 

Hence, 



dF_ _ 2r( f 
~d7 ~ 



fG {3) Mm{d - l)r~ d . (28) 



By substituting equations (27) and (28) evaluated at r = r into equation 
(26), we obtain 

d ~ 1 3 

— < 7> (29) 

or d < 4. This means that the maximum value of the spatial dimension d is 
3, and it is the upper bound we were looking for. 



4 CONCLUDING REMARKS 

We have shown that the spatial dimension must satisfy d < 4 to have stable 
circular orbits. This means that for the planetary motion in a space-time 
of five dimensions or higher, the circular orbits become instable. This result 
could be derived because the gravitational force in d dimensions possesses 
the same mathematical form as those for the central potentials studied by 
Bertrand. Although the problem we have treated in this paper is a texbook 
one [10], as we have shown above, its solution is non-trivial. Moreover, as 
was emphasized in the introduction, the concept of stability is studied in 
many current research areas of physics. Being this paper a comprehensive 
and fully detailed work, it can be a helpful pledge for the students to make 
a glance for advanced topics. 
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